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ABSTRACT 


This  report  presents  the  results  of  a  study  directed  toward  developing  a  reliability- 
prediction  technique  for  monolithic  integrated  circuits.  A  prediction  mode)  that  expresses 
reliability  as  a  function  of  device  screening,  stapling,  system  burn-in  test,  and  field  operate 
time  was  developed.  The  equation  is  based  on  data  taken  from  military,  space,  and 
commercial  application  of  integrated  circuits.  The  report  presents  the  rationale  that  led  to 
formation  of  the  equation  and  describes  its  use  and  methods  of  application. 


SUMMARY 


This  report  presents  a  reliability-prediction  equation  for  monolithic  integrated  circuits. 
The  equation  was  derived  through  multiple  linear-regression  analysis  of  field  data  on  sixteen 
system  types.  The  data  were  supplied  by  either  the  users  or  the  manufacturers  of  the 
systems. 

The  equation  presented  conforms  to  the  requirements  of  the  Weibull  distribution.  The 
shaping  parameter  (/3)  is  a  constant  in  the  equation,  while  the  scaling  parameter  (a)  is  a 
function  of  device  screening  and  sampling  tests  plus  system  bum-in  time.  The  equation 
shows  that  the  hazard  rate  continuously  decreases  with  increasing  field  operating  time, 
independently  of  the  predictor  variables.  However,  the  magnitude  of  the  hazard  rate  at  any 
specific  time  is  a  function  of  the  predictor  variables. 

The  equation  explains  the  variability  in  the  source  data  quite  well,  as  indicated  by  a 
coefficient  of  determination  (R2  )  value  of  0.87.  The  predictor  variables  are  considered  to  be 
of  the  ideal  type,  since  they  can  be  changed  (with  resultant  changes  in  reliability)  without 
restricting  the  device  or  system  designer’s  freedom  of  action. 
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CHAPTER  ONE 


INTRODUCTION 


This  report  presents  the  results  of  an  effort  directed  toward  the  development  of  a 
reliability-prediction  technique  for  monolithic  integrated  circuits;  it  reports  on  the  second 
phase  of  a  two-phase  project.  The  first  phase,  reported  in  ARINC  Research  Publication 
912-01-1-1002,  dated  October  1969,  was  concerned  with  determining  the  availability  of  IC 
reliability  data  and  the  feasibility  of  forming  a  prediction  equation  that  expressed  reliability 
as  a  function  of  design,  test,  and  application  factors.  The  initial  study  concluded  that  the 
procedure  was  feasible  and  that  data  were  available  to  support  the  development.  In  the 
current  study,  the  necessary  .data  were  acquired  and  a  prediction  technique  was  developed. 

Chapter  Two  of  this  report  describes  the  activity  that  led  to  the  formation  of  the 
prediction  model,  presents  the  prediction  model,  and  discusses  the  general  model 
characteristics.  Chapter  Three  illustrates  the  application  of  the  technique.  Chapter  Four 
provides  conclusions  resulting  from  this  study  and  recommendations  for  future  development 
activity. 
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CHAPTER  TWO 

PREDICTION-TECHNIQUE  DEVELOPMENT 


In  the  report  on  Phase  I  of  the  two-phase  study,*  we  concluded  that  sufficient  data 
were  available  to  support  the  development  of  a  reliability-prediction  technique  for 
integrated  circuits.  Further,  a  plan  for  developing  such  a  technique  was  proposed.  The  plan 
specified  da^i-collection  and  data-analysis  activities.  The  analysis  effort  was  to  depend 
heavily  on  the  use* of  multiple-regression  techniques,  which  were  expected  to  be  useful  for 
determining  the  effect  of  each  of  the  factors  influencing  device  reliability. 

To  clarify  the  nature  end  the  goals  of  regression  analysis,  a  few  preliminary  comments 
are  appropriate.  Regression  is  a  statistical  method  for  analyzing  data  to  obtain  values  of 
parameters  in  a  model  that  relates  a  dependent  variable  to  a  set  of  independent  (regressor, 
prediction)  variables.  Two  prime  applications  are  (1)  to  obtain  a  model  and  associated 
parameter  values  that  will  provide  the  best  description  or  summary  of  the  data  obtained, 
and  (2)  to  use  the  data  obtained  to  fit  the  parameters  of  a  given  general  model  for  use  in 
future  analyses  and  predictions. 

In  this  study  the  second  application  is  of  major  concern  since  we  are  interested  in 
developing  a  technique  for  predicting  the  reliability  of  integrated-circuit  devices.  As  a  result 
of  the  Phase  I  effort,  a  number  of  variables  that  are  possible  influences  on  IC  reliability  were 
identified.  Although  a  specific  model  for  relating  these  parameters  to  a  device-reliability 
measure  was  not  formulated,  two  actions  concerning  such  a  model  were  taken  on  the  basis 
of  engineering  analysis.  These  can  be  summarized  as  follows: 

(1)  It  was  hypothesized  that  the  underlying  failure  mechanisms  of  integrated -circuit 
devices  result  in  a  decreasing  hazard-rate  function  over  the  time  interval  of  general 
interest. 

(2)  The  sign  of  the  partial  derivative  of  the  reliability  measure  with  respect  to  each  of 
the  possible  prediction  parameters  was  established.  This  is  equivalent  to  establishing 
the  signs  of  the  regression  coefficients  in  the  model. 

One  immediate  restriction  concerns  the  available  data.  Suppose,  for  example,  that  the 
amount  of  acid  not  rinsed  away  prior  to  packaging  has  a  significant  effect  on  device 
reliability.  This  factor  cannot  be  directly  observed;  thus  a  related  secondary  variable  that 
can  be  measured,  such  as  device-screening  procedures  in  this  case,  must  be  used. 

Although  the  use  of  such  secondary  variables  may  clearly  weaken  the  predictive  power 
of  a  regression  model,  care  in  selecting  them  increases  the  chances  for  finding  a  satisfactory 
relationship. 


'Initial  Study  for  the  Dcvr'opmcnt  of  a  Reliability-Prediction  Technique  for  Monolithic  Integrated  Circuits, 
AR1NC  Research  Publication  912-01-1-1002,  October  1969. 
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Similarly,  since  the  data  in  this  study  were  not  obtained  through  strictly  controlled 
testing,  some  data  samples  are  lacking  one  or  more  important  variables,  making  it  necessary 
to  omit  such  variables  from  the  analysis.  The  goodness  of  fit  is  thereby  diminished,  but  if 
i  the  resulting  model  still  appears  significant  (as  measured  by  any  of  several  statistical 
measures),  it  is  possible  that  even  further  improvements  can  be  made  through  better  data 
collection. 

To  summarize,  regression  analysis  provides,  from  experimental  data,  the  coefficients  for 
the  best  fit  of  a  given  model  (e.g.,  by  the  least-squares  criterion).  It  does  not  provide  the 
optimum  model  form.  Model  selection  or,  at  least,  the  formulation  of  hypotheses  implying 
general  model  characteristics  is  the  responsibility  of  the  analyst. 

The  regression  approach  was  implemented  during  the  Phase  II  effort  described  in  this 
report.  The  following  subsections  describe  the  prediction  parameters  considered,  the  data 
sources,  the  analysis  of  the  collected  data,  the  regression  analyses,  and  analysis  results. 

2.1  PREDICTION  PARAMETERS 

The  Phase  I  effort  identified  several  variables  (listed  in  the  Phase  I  report)  that  could  be 
reasonably  expected  to  influence  reliability  significantly.  Information  on  all  of  these 
variables  was  sought  from  the  data  sources  contacted. 


2.2  DATA  SOURCES 

Industry  and  government  sources  were  contacted  to  obtain  data  samples  that  could 
provide  the  parameters  identified  in  the  Phase  I  report,  as  well  as  corresponding  field 
operating  times  and  numbers  of  failures.*  Table  1  lists  the  agencies  that  provided  data  for 
the  study;  because  of  inadequacies,  the  submittals  of  some  of  these  agencies  were  not  used. 
The  final  sample  represented  17  different  systems,  with  completely  unique  test  or 
application  situations,  containing  95  different  device  types. 


2.3  DATA  ANALYSIS 

The  collected  data  were  examined  in  the  light  of  the  parameters  of  interest,  the  data 
quantity,  and  the  variability  present  in  the  data.  Table  2  summarizes  the  specific  parameters 
of  interest  and  identifies  the  parameters  for  which  both  the  quantity  and  variability  of  data 
were  sufficient  for  analysis. 

1 

As  shown  in  Table  2,  examination  of  the  data  indicated  that  appropriate  information 
was  not  available  for  some  parameters.  For  example,  the  data  for  electrical  derating 
generally  stated  a  minimum  value  that  would  be  tolerated  as  opposed  to  the  actual  derating 
employed.  This  parameter  was  therefore  not  analyzed  further. 

For  other  parameters,  while  the  data  were  ample,  there  was  little  variation  among  the 
sources.  The  “epitaxial  technique”,  for  example,  was  use  d  in  nearly  all  cases,  so  that  the 
effect  of  not  employing  it  could  not  be  determined  from  the  available  data.  Such  parameters 
were  therefore  not  included  in  the  subsequent  analyses. 


♦For  additional  details  on  the  data  sources,  see  the  Phase  I  report,  ARINC  Research  Publication  912-01-1- 
1002,  October  1969.  This  report  is  available  at  DDC  (AD  702400). 
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Table  1.  ORGANIZATIONS  THAT  SUPPLIED  DATA 


Organization 

Location 

System 

Millions  of 
Part  Hours 

Raytheon 

Waltham,  Mass. 

Apollo 

140.8 

Reliability 

Analysis  Center 

IITRI-Chicago, 

Ill. 

Several 

5300 

Northrop 

Nortronics 

Hawthorne,  Cal. 

C5A  —  Inertial 

42.8 

Sperry 

Gyroscope 

Great  Neck,  N.Y. 

A/D  Converter 

8.7 

Loral 

Electronics 

Bronx,  N.Y. 

ECM 

3.2 

McDonnell 

Douglas 

St.  Louis,  Mo. 

Collision 

Avoidance 

1.3 

NSA 

Ft.  Meade,  Md. 

N/A 

25 

Litton 

Van  Nuys,  Cal. 

Communication 

12 

APL 

Silver  Spring,  Md. 

Satellite 

30 

IBM 

Oswego,  N.Y. 

Computers 

270 

Astronautics 

Milwaukee,  Wis. 

F-lll  -  MK  II 

1.19 

Radiation 

Melbourne,  Fla. 

AN/ASW-25A 

22.5 

Boeing 

Seattle,  Wash. 

Several 

G.E. 

Utica,  N.Y. 

3  Systems 

103.7 

UNIVAC 

St.  Paul,  Minn. 

Several 

687 

Electronic 

Communications 

Inc. 

St.  Pete,  Fla. 

10.6 

Sanders 

Nashua,  N.H. 

F-lll  Avionics 

2.8 

Litton 

Woodland  Hills, 
Calif. 

Inertial  Navigation 

120 

TRW 

Redondo  Beach, 
Calif. 

Minuteman 

10,000 

Table  2.  SUMMARY  OF  AVAILABILITY  AND  SUITABILITY  OF 
DATA  ON  PREDICTION  PARAMETERS 


Parameter 

Data  Available 
on  Most  Samples 

Data 

Variability 

Adequate 

Parameter 
Used  in 
Regression 

Design 

Basic  circuit  class 

Specific  circuit  class 

Specific  circuit  function  (complexity) 

Type  of  transistor 

Number  of  leads 

Number  of  interconnection  levels 

Buried  layers 

Guard  rings 

Epitaxy 

Isolation 

Passivation 

Number  of  diffusions 
Metalization  system 
Die  bonding  method 
Type  of  bond 
Package  type 
Device  lot  or  date  code 


Inspection  and  Test 


Device  screening 

X 

Device  sampling 

X 

Equipment/system  burn-in 

X 

Application 


Mounting  method 
Electrical  derating 
Transient  protection 
Thermal  stress 
Hermetic  equipment 
Next  assembly 

Application  for  which  data  were 
accumulated 

Type  of  system  m  which  used 
Required  equipment  reliability 
Calendar  operating  time 
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Following  selection  of  the  parameters  for  which  suitable  information  was  available, 
preparation  for  the  regression  analysis  was  initiated.  This  type  of  analysis  involves 
determining  which  values  of  an  equation’s  coefficients  best  fit  the  equation  to  the  observed 
data.  To  accomplish  this,  it  is  necessary  first  to  develop  the  form  of  the  equation  and  the 
measures  of  the  independent  variables.  In  some  cases,  Jiere  is  no  difficulty  in  quantifying 
the  independent  variables,  because  the  parameters  are  already  expressed  in  satisfactory 
terms.  For  example,  the  number  of  leads  can  be  entered  directly  as  an  independent  variable 
in  the  regression  equation.  * 

For  other  parameters,  a  binary  type  of  indication  is  appropriate.  A  “zero”  or  “one” 
indication  can  be  employed  to  show  that  an  activity  was  or  was  not  accomplished,  or  that 
one  of  two  or  more  possible  methods  was  employed. 

For  other  parameters,  however,  it  is  necessary  to  develop  a  unique  quantification 
method.  For  this  study,  screening  and  sampling  data  were  given  considerable  attention  in 
this  respect. 

Device  screening  may  consist  of  many  separate  activities  -  e.g.,  bias  bum-in, 
temperature  cycling,  and  acceleration.  If  a  large  enough  data  sample  is  available,  each  of 
these  activities  can  be  incorporated  as  a  separate  variable  in  the  regression  equation. 
However,  the  data  limitations  of  this  study  dictated  that  the  extent  to  which  the  several 
activities  were  employed  be  reflected  in  only  one  or  two  variables.  This  was  accomplished 
by  ascribing  weights  to  the  activities  to  indicate  the  estimated  relative  influence  of  each 
activity  on  reliability.  Four  engineers,  knowledgeable  in  integrated  circuits,  individually  and 
independently  assigned  relative  weights  to  the  nine  activities  considered  most  important. 
The  assigned  weight  indicated  the  individual’s  estimate  of  the  relative  importance  of  that 
activity  in  assuring  the  delivery  of  a  reliable  device. 

Each  estimator  was  permitted  to  use  any  scale  he  chose.  His  estimates  were  later 
normalized  on  a  percentage  basis  so  that  they  could  be  combined  with  the  estimates  of  the 
other  individuals.  Table  3  shows  the  normalized  weights  assigned  to  the  several  activities  by 
each  estimator  and  the  average  of  the  weights  assigned  for  each  activity. 


Table  3.  WEIGHTING  FACTORS  FOR  SCREENING  ACTIVITIES 


1 - 

Screening  Activity 

Factors  Assigned  by  Estimators 

■ 

•Estimator 

A 

Estimator 

B 

Estimator 

C 

Estimator 

D 

Operating-Life  Burn-in 

mm 

16.7 

30.7 

21.6 

22.2 

Bias  Burn-In 

■99 

11.4 

19.2 

0 

12.1 

Pre-Cap 

13.7 

16.7 

15.4 

17.6 

15.8 

Thv..mal  Shock 

10.0 

11.4 

11.5 

27.0 

15.0 

Acceleration 

11.8 

20.2 

9.6 

16.2 

14.4 

Temperature  Cycling 

10.0 

4  4 

7.7 

2.7 

6.2 

Seal 

10.0 

3.8 

10.8 

9.6 

X-ray 

5.9 

1.9 

1.4 

2.2 

Post-Cap 

2.0 

0 

2.7 

1.8 

During  the  regression  analyses,  runs  were  made  in  which  these  weights  were  varied 
considerably  to  provide  a  measure  of  weighting  sensitivity.  No  significant  differences 
resulted,  however,  and  the  original  weighting  scheme  was  retained. 

With  a  weight  assigned  to  each  activity,  it  was  then  necessary  to  provide  a  basic  raw 
score  to  which  the  weighting  factor  could  be  applied.  Generally,  a  raw  score  of  either  zero 
or  unity  was  assigned  —  unity  if  the  activity  was  performed  in  a  relatively  standard 
manner,  zero  if  the  activity  was  not  performed  or  if  it  was  performed  in  a  substandard 
manner.  For  some  activities,  e.g.,  hermetic-seal  tests,  a  value  between  zero  and  unity  was 
assigned  for  partial  performance  of  the  activity. 

For  device  bum-in,  an  expression  involving  both  test  duration  and  test  temperature  was 
needed.  The  temperature  relationship  was  developed  by  regression  analysis  of  data  supplied 
by  the  Boeing  Company;*  an  exponential  time  factor  was  then  included.  The  derivation  of 
the  complete  expression  is  given  in  Appendix  A. 

The  results  of  the  weighting  and  quantification  of  the  screening  variable  will  be  shown 
in  Section  3.1. 

To  develop  sampling-test  scores,  the  same  approach  employed  for  screening  was  applied; 
i.e.,  the  same  raw  scoring  and  weighting  factors  were  used.  Two  additional  activities 
(high-teinperature  storage  and  life  test)  were  included  for  the  sampling  score. 

When  the  effects  of  sampling  tests  were  considered  initially,  it  was  suggested  that  the 
specified  lot-acceptancc  sampling  plan  would  influence  the  resulting  reliability.  The  measure 
employed  to  reflect  the  sampling  plan  -as  the  Lot  Tolerance  Percent  Defective  (LTPD),  a 
standard  quality-control  statistic.  In  the  regression  analysis  that  followed,  several  methods 
of  including  the  LTPD  were  tried.  The  most  satisfactory  results  were  obtained  when  the 
weighted  score  was  divided  by  the  LTPD.  The  sampling-test  scoring  results  are  given  in 
Section  3.1. 

2.4  REGRESSION  ANALYSES 

The  approach  used  in  this  study  was  to  begin  the  regression-analysis  work  with  a  large 
group  of  variables  and  use  the  results  of  the  analyses  to  guide  successive  selection  of 
variables  and  model  forms.  Preference  was  always  given  to  the  more  primary  variables,  and 
the  model  forms  used  were  influenced  by  the  hypothesis  that  the  device  failure  rate 
decreased  with  time. 

During  the  course  of  the  program,  95  separate  runs  were  made  before  the  results  were 
considered  satisfactory.  While  each  run  cannot  be  discussed  in  detail,  the  analyses  are 
summarized  in  Table  4,  with  the  runs  grouped  in  18  sequences. 

The  table  includes  the  variables  that  were  found  to  contribute  significantly  to  the 
predictive  ability  of  the  equation.  It  should  be  noted  that  all  of  the  variables  listed  in  Table 
2  for  which  suitable  data  were  available  were  used  in  the  initial  ru  is»  The  “Significant  Input 
Variables”  column  in  Table  4  lists  tho.-j  variable.'  which  were  not  rejected  by  the  analysis. 


*Reliability  Chaiactcriration  and  Prediction  of  Integrated  Circuits,  RADC  TR-70-232,  prepared  by  the  Boe¬ 
ing  Company,  Aerospace  Group,  November  1970. 
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Table  4  S  I'M  MARY  OK  REGRESSION  ANALYSIS 


rrun-ifurnutiun  tncd  on 
Dependent  Variable 


6  I  (Same  as 


Weighting 
Functions  Used 


Device  Family 

Packaging 

Complexity 

Reliability  Requirem  -tt 
Sampling  Scon- 
Screen  mg  Score 
Application 

(Same  as  1  except  for 
Pac  Waging) 


Device  Burr-In 
Sampling  Score 
Screening  Score- 
System  Test 


(Same  as  4) 


Equal  0.28-0.70  j  Eliminated  pa  Paging  due  to 

VES~5tT»KM«  I  n“ml" ' ,,r  T0  5 

Bn  (operating 
time) 


Bn  (operating  0.31-0  .6$^.  Lumw\aV.i  n.vrpvs  with  no 
bme)  recorded  failure* *  Reduo-d 

complexity  to  two  categoric' 


7 

(Same  as  -1) 

8 

8 

System  Test  Houis 
Screening  Score 

Sampling  Score 

Average  Calendar  Time 
in  Field 

* 

9 

(Same  as  8) 

15 

Bn  (operating 
time) 


Bn  (operating 
time) 


Cn  (operating 
bme) 


0.44-0.68  Removed  device  bum  m  fron 
screening  an  J  added  ty»b-n 
test 


Ch;uv  .  in  input  variable- 


Assumed  one  failure-  for  sam¬ 
ples  with  no  recorded  failures 


Bn  (operating  0.53-0.74 
time) 


Bn  (System  Test  Hours) 
Bn  (Average  Calendar 
Time  in  f  ield) 


Bn  (Screening) 
Bn  (Sampling) 


Bn  (Sampling  ♦  Screen¬ 
ing) 


Bn 

Equal 

No.  of  Failures 

s/NoTofT  ai  lures 
Bn  (operating 
hour>) 

0.51  0.95 

Bn 

Equal  . 

No  of  failures 
v^No.  of  f  ailures 
Bn  (operating 
hours) 

0.42-0.96 

Used  only  samples  with  > 

5  X  10*  oj>orJtinK  hours 
Tried  different  weighting 
•chcmi-s  for  »•  reeling  «'or«- 


Major  system  totals.  V.‘H 
lower  50'-  confidence  limit 
for  tero  failure  sample* 


Equal 

0.53-0  82 

Add  one  t* 

Mo.  of  Failures. 
y/So.  of  Failures 
(n  (operating 
hour--) 

Operating  hours 

of  failure- 

Screening  Score 
Sampling  Score 
System  Test  Hours 


System  Tret  Hour? 
Screening  Score 
Sampling  Score 
Average  Calendar  Time 
m  Field 


0.50-0.75 

Eliminated  samples  with 

xero  failures 

umNt  of  failures 


No.  of  Failures  0.65-0.79  Eliminate-)  sa: 

Bn  (operating  failures 

time) 

Equal 

y/So.  of  Failure* 

In  (op.  tune)  X 
No.  of  Failures 

Equal  0.46-0.94  Lower  50 ‘  confidence  limit 

UKsrsmsxi  ,or  “"■i-1"  *,ih  '-■ro  f*-:urp> 

Bn  (op  time)  X 
No.  of  Fa;  Wes  ' 

No  of  Failures  j  0.67-0  S3  j  Lower  50 

jssrarasi;  <«’  » 


I  0  67-0  S3  Lower  50  .  .Mif.Ji-ncc  L-iiit 

!  for  samp!-  s  with  *cro  failures 

|  New  sampling  s<  <»n» 

I  0.72-0  77  I  Lower  50~i  ccnfidi  m  e  limit 


‘Natural  logarithm . 


The  tabulation  also  gives  the  number  of  runs  made  in  each  group.  This  number  is  often 
indicative  of  the  transformations  made  in  the  variables  and  the  relative  weights  accorded 
each  data  set.  In  the  matter  of  transformations,  it  is  sometimes  useful  to  vary  the  form  of 
the  basic  equation,  for  example,  from  a  sum  to  a  product  function.  This  is  easily 
accomplished  by  using  logarithms  of  the  variables. 

It  is  often  advisable  to  './eight  the  data  sets  in  terms  of  the  size  of  the  data  base  so  that  a 
relatively  small  sample  will  not  have  a  disproportionate  effect  on  the  overall  result.  Since  the 
amount  of  information  in  a  test  sample  is  related  to  the  number  of  failures,  this 
characteristic  is  often  used  in  weighting  activities. 

The  final  statistic  noted  in  Table  4  is  the  range  of  R2  values  experienced  in  the  group  of 
runs.  The  R2  (coefficient  of  determination)  value  reflects  the  amount  of  variation  in  the 
data  that  is  explained  by  the  developed  equation.  It  is,  perhaps,  the  best  indicator  of  how 
well  the  equation  fits  the  data  being  analyzed. 

% 

Table  4  also  includes  specific  remarks  on  significant  .or  peculiar  characteristics  of 
regression  runs  in  each  group. 

In  the  initial  regression  runs,  many  input  parameters  were  considered,  and  those  which 
appeared  to  have  the  greatest  effect  on  device  reliability  were  identified.  The  number  of 
observations  varied  because,  of  voids  in  the  data  set.  For  example,  54  samples  contained  a 
description  of  the  package  type.  When  it  was  observed  that  the  failure  rate  did  not  depend 
greatly  on  the  package  type,  the  packaging  parameter  was  omitted  from  consideration,  thus 
permitting  an  additional  2G  samples  to  be  considered.  This  procedure  continued  until  the 
final-screening,  sampling,  burn-in,  and  field  operate  parameters  were  identified  as  those 
having  the  greatest  impact  on  the  observed  failme  rates. 

The  final  set  of  runs  (sequence  numbers  16  through  18)  led  directly  to  the  results 
presented  in  this  report.  For  trese  runs,  the  samples  for  which  no  failures  were  reported 
were  assigned  0.69  failure  (the  0.69  is  calculated  from  the  lower  50-percent  confidence  limit 
of  the  chi-squared  distribution*).  In  addition,  the  reciprocal  of  the  LTPD  was  introduced 
for  the  calculation  of  the  sampling  score.  The  data  points  were  weighted  according  to  the 
number  of  failures  and  the  square  root  of  the  number  of  failures,  and  two  basic  forms  for 
the  regression  equation  were  assumed.  All  of  the  later  regression  runs  indicated  a  good  fit  to 
the  data  (i.e.,  high  R2 );  it  was  therefore  possible  to  introduce  engineering  judgment  in 
selecting  -the  final  regression  formulation.  It  was  decided  that  the  form  of  regression 
equation  that  matched  the  Wcibull  distribution  would  have  the  greatest  general  applica¬ 
bility.  It  was  also  decided  that  the  data  samples  should  be  weighted  according  to  the  square 
root  of  the  number  of  failures  so  that  the  observations  with  higher  numbers  of  failures 
would  have  a  larger  but  not  overwhelming  effect  on  the  regression  equation. 


♦In  a  strict  sense,  this  approximation  would  be  applicable  only  if  the  device  failure  rate  were  constant. 
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2.5  DISCUSSION  OF  RESULTS 


2.5.1  Regression  Equation 

The  equation  resulting  from  the  regression  analysis  that  was  considered  most 
appropriate  and  useful  for  prediction  purposes  is* 

4*  =  76  877 1/3  e0.025  Sc  +  0.00095  Sa  +  0.0064  L 

where  . 

4>  =  total  field  operating  time  divided  by  total  failures 
tj  “  system  burn-in  time  (hours) 
t2  "  field  operate  time  (hours) 

Sa  =  acceptance  sampling  score 
Sc  m  screening  score 


This  equation,  based  upon  17  observations,  exhibits  a  coefficient  of  determination,  R* , 
of  0.87. 

The  observations  represented  a  combination  of  the  95  device  types  into  17  groups  for 
each  of  which  the  parameters  of  the  equation  were  the  same  —  that  is,  all  devices  in  each 
of  the  17  groups  achieved  the  same  screening  and  sampling  scores  and  experienced  the  same 
burn-in  time.  (The  devices  in  each  group  achieved  the  same  scores,  of  course,  because  they 
were  used  in  the  same  system.)  Combination  in  this  manner  also  minimized  the  statistical 
difficulty  encountered  when  no  failures  were  observed  in  a  particular  data  set.**  The 
grouping  of  data  in  this  manner  resulted  in  a  zero-failure  experience  for  only  one  of  the  17  ' 
groups. 

A  preliminary  evaluation  of  the  predictive  ability  of  the  equation  was  performed.  The 
equation  was  used  to  predict  4>  for  the  devices-that  provided  the  data  from  which  the 
equation  was  developed.  The  predicted  4>  for  each  of  the  data  samples  was  then  plotted 
against  the  observed  4>  for  that  sample.  The  results  are  shown  as  a  scattergram  in  Figure  1. 

While  this  can  hardly  be  considered  a  rigorous  test  of  the  equation,  it  is  encouraging  to 
note  that  with  the  exceptions  of  the  two  points  that  lie  well  abrve  the  45°  line,  which 
indicates  perfect  agreement,  there  is  relatively  little  deviation  from  the  line.  Further,  one  of 
the  two  points  with  large  deviations  represents  a  sample  that  had  no  failures,  and  the  other 
point  represents  a  sample  with  only  one  failure.  In  the  analysis,  they  carried  little  weight 
since  the  equation  was  derived  with  the  samples  weighted  according  to  the  square  root  of 
the  number  of  failures  experienced.  (A  value  of  0.69  failure  was  assumed  for  the  zero-failure 


♦The  exponent  of  t2  war  originally  calculated  to  be  0.3GC2.  To  facilitate  the  use  of  the  prediction 
equation,  a  ridge  regression  analysis  was  employed  to  adjust  the  exponent  to  1/3.  The  adjustment  had  an 
insignificant  effect  or;  the  R5  value,  decreasing  it  from  0.8736  to  0.8735.  The  ridge  regression  analysis  is 
discussed  in  Hoerl,  A.E.  and  Kennard,  II. W.,  “Ridge  Regression:  Biased  Estimated  for  Nonorthogonal 
Problems”,  Tcclmotnc  tries,  Vol.  12,  No.  1,  February  1970,  pp  55-68. 

**It  will  be  recalled  that  many  of  the  “remarks”  on  Table  4  referred  to  various  methods  for  treating  “zero 
failure”  data  sets. 
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Actual  <!>  (Hours  +  Failures) 


sample.)  Therefore,  those  two  points  had  less  influence  on  the  equation  than  suggested  in 
Figure  1,  as  evidenced  by  the  high  value  of  R1  associated  with  the  equation. 

The  prediction  error  is  less  than  2  to  1  for  more  than  80-percent  of  the  samples  and  less 
than  1.5  to  1  for  the  majority  of  the  samples;  the  dashed  lines  in  Figure  1  represent  these 
error  ratios. 

2.5.2  Derivation  of  Reliability  Function 

The  model  indicates  that  the  reliability  of  integrated  circuits  is  a  function  of  field 
operating  time.  The  prediction  equation  for  4>  can  be  used  to  obtain  the  parameters  of  a 
Weibull  distribution  and,  consequently,  a  reliability-prediction  equation.  Such  a  derivation  is 
described  in  Appendix  B.  The  reliability  function  derived  is  presented  below,  together  with 
the  range  of  parameter  values  supported  by  the  data: 


where 

Rw(t»)  is  read  as:  Rw  is  a  function  of  tj 
K  =  76  877  Sc  +  0.00095  Sa  +  0.0064  t] 

Sc  "  screening  score  0  <  Sc  <  70 
Sa  -  sampling  score  0  <  Sa  <  2438 
tt  «  system  bum-in  time  0  <  tj  <  311 
tj  -  field  operate  time  96  <  t2  <  14,600 

Three  particularly  significant  points  can  be  made  about  the  equation: 

(1)  The  equation  reveals  a  decreasing  failure  rate.  This  has  long  been  suspected  by 
many  in  the  field,  but  previous  prediction  techniques  have  not  accounted  for  this 
factor. 

(2)  With  the  exception  of  field  operating  time,  all  variables  in  the  equation  can  be 
controlled  within  wide  limits  during  production.  Thus  if  the  predicted  reliability  is 
not  acceptable,  the  equation  implies  that  changes  can  be  made  to  improve  the 
reliability.  Further,  since  these  changes  would  be  concerned  only  with  testing, 
neither  the  device  nor  the  system  designer’s  freedom  of  action  would  be  affected. 

(3)  The  equation  implies  that  the  underlying  failure  mechanisms  are  such  that  the 
Weibull  shape  parameter  is  not  influenced  by  the  predictor  variables  but  the  scale 
parameter  (a)  is.  For  the  simple  exponential  distribution,  for  example,  this  means 
that  the  constancy  of  the  hazard-rate  function  is  independent  of  the  predictors  but 
the  mean  life  is  not. 

Thus  ARINC  Research  believes  that  the  equation  presented  is  the  most  accurate  and 
useful  one  that  has  been  developed  for  integrated  circuits. 


13 


CHAPTER  THREE 


TECHNIQUE  APPLICATION 


This  chapter  presents  the  processes  to  be  employed  for  determining  the  equation 
parameters,  gives  an  example  of  the  application  of  the  equation,  and  provides  comments  on 
the  areas  of  applicability  and  the  limitations  of  the  equation. 

3.1  DETERMINING  THE  EQUATION  PARAMETERS 

The  first  and  most  difficult  step  in  using  the  prediction  equation  is  to  determine  the 
equation  parameters.  The  difficulty  lies  in  the  device  sampling  and  screening  scores;  system 
bum-in  time  and  field  operating  time  are  used  directly  in  the  equation. 

Table  5  illustrates  the  derivation  of  the  device  sampling  and  screening  scores.  A  raw 
score  is  derived  for  each  test  on  the  basis  of  the  specific  test  conditions  and  requirements. 
For  screening,  each  of  these  raw  scores  is  then  multiplied  by  the  appropriate  weighting 
factor.  (The  derivation  of  these  weighting  factors  was  described  in  Section  2.3.)  The  sum  of 
these  products  is  the  desired  screening  score. 

For  the  sampling  score,  an  additional  step  is  required:  the  product  of  the  raw  score  and 
weighting  factor  is  divided  by  the  Lot  Tolerance  Percent  Defective  (LTPD).  The  individual 
scores  for  the  tests  are  then  added  to  obtain  the  final  sampling  score. 

The  specific  scores  shown  in  Table  5  correspond  to  MIL-STD-883,  Level  B.  Since  this 
standard  is  widely  employed,  the  sampling  and  screening  scores  that  would  be  used  if  Levels 
A,  B,  or  C  of  the  standard  were  employed  are  presented  as  follows: 


MIL-STD-883 

Screening 

Sampling 

Levels 

Score 

Score 

A 

102.5 

960 

B 

67 

713 

C 

40 

368 

While  guidance  on  scoring  the  screening  and  sampling  variables  is  provided  by  Table  5,  it 
is  obvious  that  not  all  possibilities  can  be  covered  specifically.  If  the  particular  test  being 
scored  is  not  one  of  the  standards  and  the  conditions  or  criteria  are  not  covered  by  Table  5, 
engineering  judgment  will  have  to  be  used.  Some  comment  on  two  scoring  problems 
encountered  in  this  study  and  how  they  were  resolved  may  assist  the  engineer  in  ranking  his 
judgments. 


Preceding  page  blank 
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Table  5.  SCORE  SHEET  FOR  SCREENING  AND  SAMPLING  VARIABLES 


Test 

Screening  Score 

Sampling  Score 

Raw  ,  ,  Weighted 

Score  x  V/ei2ht  =  Score 

Raw  „  Weighted 

Score  *  V*  eight  LTPD  =  Score 

1.  Pre-Cap  Visual 

•  For  M1L-STD-883  Level  Acr  11  or  equivalent, 
let  RS  =  1.0 

•  For  M1L-STD-883  Level  C  or  equivalent,  let 
RS  =  0.75 

•  For  neither  of  the  above,  let  RS  =  0 

1.0  X  16  =  16 

OX  16  v  -  =  0 

2.  Post  Cap  Visual 

See  #1  for  raw  scores 

10X2  =  2 

1.0  X  2  v  0.15  =  13.33 

3.  Thermal  Shock 

•  From  at  least  -55  to  +125<>C,  let  RS  =  1.0 

•  For  no  test,  let  RS  =  0 

OX  15  =  0 

1.0  X  15  4  0.15=  100 

4.  Acceleration 

•  For  gold  leads,  20000G’s,  Y,  and  Y2  direc¬ 
tions,  let  RS  =  1.0 

•  For  gold  leads,  20000G’s,  Y|  only,  let  RS  = 
0.7 

•  For  aluminum  leads,  30000G's,  Y,  and  Yj , 
let  RS  -  1.0 

-  For  aluminum  leads,  30000G's,  Y,  only,  let 
RS  -  0.7 

-  For  none  of  the  above,  let  RS  =  0. 

0.7  X  14  =  9.8 

0.7  X  14  v  0.15  =  65.3 

5.  Temperature  Cycling 

•  From  at  least  -  55  to  125°C,  let  RS  =  1.0 

•  For  no  lest,  let  RS  =  0 

1.0  X  6  =  6 

6.  Seal 

•  For  MIL-STD  fine  and  gross,  let  RS  =  1.0 

■  For  fine,  only,  let  RS  =  0.3 

■  For  no  test,  let  RS  =  0 

1.0  X  30  =  10 

1.0  X  10  4  0.15  »  66.7 

7.  X-ray 

•  For  gold  leads,  let  RS  -  1.0 

■  For  aluminum  leads,  let  RS  =  0.3 

•  For  no  test,  let  RS  =  0 

0X2  =  0 

OX  2r-*=0 

8.  Burn-in  (Opera!  ional);  Compute: 

•  RS  =  104.6(l-o't/56)  •  (c'3-96  +  002T) 

•  where  t  =  Time  (hours) 

T  *  Temperature  (®C) 

23  X  1.0  =  23 

OX  1.0  4  -  =  o 

9.  Bum-in  (Bias);  Compute  RS  ns  in  #8 

OX  0.5  =  0 

31.6  X  0.5  4  0.10  =  158 

10.  High-Temprrature  Storage 

•  For  temperatures  of  at  least  125°C,  and  time; 

>500  hours,  let  RS  =  1.0 
■  For  no  test,  let  RS  =  0 

1.0  X  6  4  0.15  =  40 

11.  Life  Test*;  Compute: 

•  RS  =  104.6  (1-e"1/333)  (c'3  9fi  +  002T) 

•  where  t  =  Time  (hours) 

T  =  Temperature  (’C) 

23  X  1.0  v  0.10  =  230 

Combined  Screening 

Score . 66.8 

Combined  Sampling 

Score . 713.3 

Note:  RS  *■  Raw  Score. 

RS  to  be  determined  separately  for  screening  and  sampling. 


•Credit  cannot  ho  fiven  for  both  test  11  and  ted  8. 


/&> 


In  one  case,  the  device  specifications  did  not  include  any  lot  s  .npling  requirements. 
However,  for  all  screening  tests,  the  specification  required  that  a  lot  be  rejected  if  the 
percentage  failing  the  screening  test  exceeded  a  certain  level.  The  sampling  score  was  based 
on  an  LTPD  equal  to  the  percent  defective  that  would  result  in  lot  rejection.  While  this 
score  was  admittedly  somewhat  low,  we  were  unable  to  Llerive  the  exact  equivalent  LTPD, 
since  the  sample  constituted  100  percent  of  the  lot.  ' 

In  another  case,  temperature  cycling  was  specified,  but  the  range  was  -50  to  +125° C. 
Full  credit  was  given  for  this  test,  even  though  the  lower  temperature  was  5°C  higher  than 
the  minimum  specified  in  Table  5. 


3.2  THE  PREDICTION  PROCESS 


Having  developed  the  scores  for  screening,  sampling,  and  system  bum-in,  we  are  ready 
to  proceed  with  our  prediction.  For  purposes  of  illustration,  suppose  we  wish  to  predict  the 
mission  reliability  of  a  given  system  for  the  following  conditions: 

Mission  length  (£ )  =  50  hours 

Number  of  integrated  circuits  in  system  (N)  =  5,000 

Field  time  on  system  at  start  of  mission  (t* )  =  3,000  hours 

Device  screening  score  (Sc)  =  50 

Device  sampling  score  (Sa)  =  500 

System  burn-in  time  (t, )  =  200  hours 

The  basic  equation  for  calculating  the  predicted  reliability  of  integrated  circuits, 
developed  in  Chapter  Two,  is  repeated: 

_  *  2/3  #rr 
Rw{t2)  -  c  /K 

where 

Rw(lj )  is  read  as:  Rw  is  a  function  of  t2 
K  =  76  877  e+0-025  +  0.00095  Sa  +  0.0064  t] 

Sc  *»  screening  score,  0  <  Sc  <  70 
Sa  -  sampling  score,  0  <  Sa  <  2,438 
t,  »  system  burn-in  time,  0  <  t,  <  311  hours 
t2  «  field  operate  time,  96  <  t2  <  14,600  hours 

For  an  integrated  circuit  operated  from  time  zero  to  time  t2 ,  the  probability  of  no 
failure  is  given  by 


^mission  “  IWUz) 
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If,  however,  the  mission  begins  at  time  t2  and  continues  until  time  t2  +  2,  the  mission 
reliability  (assuming  an  operating  system  at  t2 )  is  well  approximated  by 


^mission 


nw(t2  +  2) 

Rw(t2) 


If  a  system  contains  N  integrated  circuits  in  series  for  reliability  purposes,  all  of  which 
were  subjected  to  the  same  screening,  sampling,  bum-in  and  field  operation  *,  the  system’s 
mission  reliability  is  expressed  by  either** 


remission  “  (Rv/(f2  ))^  or 


R 


mission 


Ry;(t2  *  2) 
Rw(t2 ) 


N 


For  our  example,  then,  the  second  of  these  equations  is  applicable 


R 


(t,  +  2) 


2/3~?N 


K 


(t2)2/3 

K 


(3050) 


2/3 


K 


(3000)2/3 

K 


5,000 


where 

K 


76  877  c^-025(50)  +  0.00095(500)  +  0.0034(200)  «  li552i000 


or 


R 


c- 0.0001355 
e-0.0001340 


e- 0.0075  .  o.9925 


♦If  integrated  circuits  have  been  subjected  to  different  procedures  or  operation,  the  product  of  the 
reliabilities  determined  fer  each  circuit  wit!  yield  the  rys'em  reliability. 

♦♦These  expressions  imply  that  the  system  is  composed  only  of  integrated  circuits.  Predictions  of  the 
reliability  of  other  components  must  be  nituV'  by  conventional  m  r.ns  and  combined  with  the  results 
obtained  here. 


i 
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To  demonstrate  the  dependence  of  reliability  on  field  operating  time,  consider  the 
reliability  of  the  same  system  for  the  same  50-hour  mission  beginning  at  t2  =  500  hours: 


R 


0/-n  5,000 
(550)2''3  5 


K 


(500) 


,2/3 


K 


0.9S68 


While  both  reliabilities  are  high,  the  mission  beginning  at  3,000  hours  has  a  higher 
probability  of  success  because  of  the  decreasing-failuie-rate  characteristic. 

The  equation  may  also  be  employed  to  estimate  the  expected  number  of  failures  over 
some  stated  period  of  operating  time.  Such  a  result  would  be  useful,  for  instance,  for 
sparcs-provisioning  purposes.  The  expected  number  of  failures  (including  replacement 
failures)  during  the  period  from  t2  to  t2  +  £  can  be  approximated  by 


E 


N 

K 


(t2  +  fi)2/3  - 


This  equation  is  based  on  approximating  the  mean  number  of  renewals  by  the  negative  of 
the  logarithm  of  the  reliability  function.  Such  an  approximation  is  quite  accurate  when  the 
reliability  is  close  to  1.0. 

3.3  TECHNIQUE  APPLICABILITY  AND  LIMITATIONS 

Special  emphasis  must  be  given  to  the  usual  cautions  associated  with  the  use  of 
empirical  equations  for  conditions  not  covered  by  the  data  from  which  they  were  derived. 
While  the  range  of  individual  parameter  values  covered  by  the  data  is  quite  wide,  certain 
critical  combinations  of  parameter  values  did  not  occur  in  the  data.  For  example,  there  were 
no  cases  in  which  both  sampling  and  screening  scores  were  zero.  Therefore,  the  equations 
presented  should  not  be  used  if  both  sampling  and  screening  scores  are  zero. 

Chapter  Two  described  the  range  for  the  variables  used  in  the  prediction  equation. 
However,  there  arc  other  items  that  we  must  logically  be  concerned  with,  such  as  device 
complexity.  The  most  complex  circuit  represented  in  the  data  was  an  eight-bit  shift  register; 
no  other  medium-  or  large-scale  integration  was  included. 

Table  G  summarizes  the  scope  of  all  characteristics  of  dat  a  that  were  used  in  developing 
the  prediction  equation.  The  reader  is  warned  that  applying  the  equation  for  conditions  not 
covered  in  the  table  is  extremely  risky.  If  the  equation  is  ur^ci  for  cases  that  are  within  the 
conditions  covered  by  the  samples  used  in  the  equation’s  development,  we  believe  that  it 
will  yield  the  most  accurate  results  possible  with  any  published  technique. 
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Table  6.  RANGE  0?  DATA  CHARACTERISTICS  SCHEDULED  IN 

DATA  SAMPLE 

Characteristic 

Barge 

Technology 

Monolithic  bipolar  circuits  only 

Package  Types 

Metal  F/P,  commie  F/P 

Number  of  Device  Leads 

10  to  14 

Circuit  Types 

Digital  and  liner:/ 

Logic  Types 

DTL,  TTL,  RTL 

Bond  Types 

T.C.  Ball,  T.C.  Y.'adje,  Ultrasonic 

Metallization-Lead  Mate  /irks 

Au-Au,  Al-Al,  Au-/.l 

Complexity 

Single  Gate  to  I'-'-'i  .Shift  Register 

Date  of  Device  Mamikeiu  v 

1964-1970 

Passivation  Types 

SjO, ,  Glass 

Application  Environment 

Missile,  Aircraft,  Spacecraft,  Ground,  Sub¬ 
surface,  Shipboard 

Qualification  Class 

Hi-Rel.,  Miliir.y,  Commercial 

Screening  Score 

0  to  70 

Sampling  Score 

0  to  2,438 

System  Burn-in  Time 

0  to  311  hours 

Field  Operate  Time 

9G  to  14,600  hours 

CHAPTER  FOUR 


CONCLUSIONS  AND  UECOIdMIiNDATlONS 


4.1  CONCLUSIONS 

The  following  conclusions  resulted  from  this  study: 

•  The  equation  developed  in  the  study  reveals  that  the  most  important  predictors  of 
reliability  of  integrated  circuits  r.  e  device-screening  and  -sampling  requirements, 
system  burn-in  time,  and  field  operate  time. 

■  The  equation  fits  the  data  from  vhieh  it  vras  derived  extremely  well,  as  evidenced  by 
the  R2  value  of  0.87  and  the  fact  that  the  majority  of  the  data  points  fall  within  50 
percent  of  the  predicted  value. 

’  All  parameters  contained  in  the  prediction  equation  are  independent  of  the  design 
process.  Therefore,  the  reliability  of  a  system  can  be  improved  by  requiring  more 
stringent  testing  (c.g.,  lower  LTPD,  longer  burn-in)  without  intr  rfeiing  with  the  basic- 
device  or  equipment  designs.  While  this  does  not  moan  that  vc’v  biliiy  is  insensitive  to 
design,  the  customer  or  manufacturer  need  not  always  resort  to  redesign  to  improve 
reliability. 

•  The  results  of  this  study  strongly  support  the  hypothesis  that  integrated  circuit 
display  a  decreasing  failure  r<  to  with  increasing  operating  tv  re  over  the  period  of 
genera1  interest.  The  study  also  indicates  that  this  process  can  be  accmalely 
represented  by  the  Weibull  distribution. 

4.2  RECOMMENDATIONS 

The  following  recommendations  are  offered: 

'  The  prediction  equation  should  be  tested  and  validated  with  an  entirely  new  set  of 
data. 

•  The  range  of  variables  and  conditions  covered  by  the  equation  should  be  expanded, 
with  emphasis  on  medium-  and  large-scale-integration  circuits. 

•  Sufficient  additional  data  should  be  obtained  to  permit  derivation  of  the  sampling 
and  sercihing  scoring  factors  from  data  rather  than  amigummil  of  the  factors  by 
judgment. 

-  Con-’-1  •  ■bon  should  Ire  given  to  thriving  similar  prediction  equations  for  hybrid 
circuit:  an!  for  discrete  diodes  ;m:I  i  mi:  Hors.  The  by!-  '  1  t  A  nil  is  consult  o  .1  »’  e 

more  in.pmtnnt  at  this  time. 

The  ncco.-  ary  data  to  support  the  actions  recommended  above  could  hr  obtained  (on 
a  single,  c— v-ful!y  designed  dela-cotl 'c' ion  j  ogram.  It.  is  rccoi.v.  :  -.r’l-d  that  such  a  prog  .mi 
be  designed  re  d  implemented  bmi.'dhlely. 


APPENDIX  A 


Tiiv3E-Tr,:.::»i:RATUUti  RobATio;  ,:r.  jibs 


The  efficiency  of  cloviv.--  'urn-in  tests  is  a  function  of  both  time  and  tempers  Hire.  It  is 
usually  considered  to  cl.r  ,  ■  much  more  rapidly  with  temperature  than  with  tin... 
Therefore,  it  is  dcsivabV  to  i'  vdoji  an  expression  that  r  -cor;- is  for  both  of  these  winblos 
and  is  more  strongly  influent-,  d  by  temperature. 

The  effect  of  ten:;  ...  fere  ’.ms  obtained  by  perform/;-."  a  >  -gresslon  analysis  on  dale, 
developed  by  Boeing.-  T.-l  le  /.-]  presents  the  dale  v. .  a  i.  tV-  analysis.  The  developed 
equation  was  as  follow.--; 

-3.96  +  0.02T 
> .  -  c  < 


where 

X  »  failure  rate  vjjeriv:d  ;•  v  1.099  hours) 
Tc  -  ambient  test  .  in.  i.t  °C 


'the  time  function  \ .„  •  >  ’  d  more  arbitrarily.  ]>•»; .  l.-nco  with  all  types  of  silicon 
semiconductor  devices  J,.  ;  d o  that  168  horns  is  ;o  ;  quef-;  burn-in  period  at  h*;.h 
temperatures.  Therefor.',  rr.  <  p  ession  was  fonn.  d  that  rc’.ii  ■  95-pcrccnt  of  its  effr  cl  in 
168  hours.  This  expression  1<  .  !  the  following  form; 


l-c-‘/56 

where  t  =  the  burn-in  time  in  her.'. 


The  time  and  temp;  hue  c  n-essions  were  combined  i  •  (?,*•  following  form; 


Wei;  lit 


K,  (1  -  e-^)  {c-a.M  +  0.02T, 


The  score  is  more  fd.nug..'  inclueneed  by  a  decrerneg  T  than  by  decreasing  1.  For 
example,  if  T  dec; a..;  f;  If- n°C  to  25° C  (a  factor  of  f>),  the  T  function  drops  by  a 

I;  -foi  of  7.4.  SimlL.  -i.  .  i"  i  d  .  vs  from  162  to  88.8  (•  ’  a  factor  cd  8),  the  t  fur.,-:-o.; 
drops  only  by  e.  factor  of  :  loci  2.1 . 

To  provide  the  fin;  1  score,  K ,  was  set  at  10  :.C,  to  y;  1  ’  a  Intel  weight  of  23  for  a  t  . ; 
i*i iic  of  168  boms  r.l  r  I  ,  .  lu.v  of  1  ?5nCh 

iW  ‘a  of  Irtcy.u  !•■(!  ( 7  Y ; ’-’i  0-?.g?, )  ■-  riv  >c,l  by  0;.  1 

1  .r,  1970. 


Company,  Aero -pace  C  o::;-.  K.  a. 


Tabic  A—l.  SAMPLE  FAILURE  HATES 
(Percent  for  1000  Kov-s)* 


l’cmpciutmc 


25°C 


50°  C 


75°C 


100°C 


125°C 


Failure  Rates 


7.0,  1,8,  1.7,  1.3,  0.75,  0.70,  0.00,  0.48,  0.35,  0.30,  0.25,  0.19,  0.12, 
0.090,  0.095,  0.030,  0.075,  0.030,  0.050,  0.04b,  0.012,  0.033,  0.030, 
0.030,  0.025,  0.025,  0.025,  0.025,  0.023,  0.021,  0.020,  0.020,  0.020, 
0.018,  0.018,  0.017,  0.016,  0.013,  0.010,  0.0090,  0.0090,  0.0080, 
0.0080,  0.0075,  0.0070,  0.0070,  0.0050,  0.0010,  0.0035,  0.0030, 
0.0028,  0.0025,  0.0020 


3.0,  2.5,  2.2,  0.65,  0.60,  0.22,  0.2.2,  0.15,  C.10,  0.08,  0.07,  0.065, 
0.0G0,  0.055,  0.052,  0.052,  0.050,  0.010,  0.025,  0.020,  0.018,  0.018, 
0.015,  0.013,  0.012,  0.011,  0.009,  0.007,  0.005,  0.0015,  0.0041, 
0.0038,  0.0024,  0.0018 

3.8,  1.1 
4.5 


22,  20,  12,  9,  8,  7,  6.5,  6.0,  4.5,  4.0,  3.8,  3.5,  2.8,  2.7,  2.7,  2.4,  2.2, 
2.2,  2.2,  1.8,  1.7,  1.7,  1.7,  1.3,  1.2,  1.0,  1.0,  1.0,  0.8,  0.7,  0.7,  0.65, 
0.56,  0.65,  0.65,  0.6,  0.6,  0.6,  0.6,  0.6,  0.5,  0.5,  0.4,  0.4,  0.4,  0.4,  0.4, 
0.35,  0.3,  0.28.  0.25,  0.25,  0.25,  0.22,  0.20,  0.18,  0.18,  0.1S,  0.15, 
0.15,  0.15,  0.14,  0.13,  0.13,  0.13,  0.12,  0.12,  0.12,  0.11,  0.10,  0.10, 
0.10,  0.09,  0.09,  0.09,  0.0S,  0.08,  0.075,  0.075,  0.075,  0.070,  0.060, 
0.060,  0.060,  0.055,  0.050,  0.050,  0.040,  0.040,  0.03S,  0.038,  0.03S, 
0.030,  0.030,  0.030,  0.023,  0.028,  0.028,  0.025,  0.022,  0.022,  0.022, 
0.020,  0.020,  0.020,  0.01S,  0.017,  0.014,  0.012,  0.011 , 0.011 , 0.010, 
0.010,  0.005,  0.005,  0.0015 


*Dala  Source:  Doein;;  Study. 


APPENDIX  D 


DERIVATION  OF  RELIABILITY  FUNCTION 


The  prediction  equation  presented  in  Chapter  Two  can  be  written  as 
Total  field  operate  time  ,  ,o 

4>  -  —  —  —  -1 - -  Kt,1/3  (1) 

Number  of  failures 


where 


l2  “ 

K  ■“ 


field  operate  time 

■*  K,Sc  +  I'm  l  j 


To  derive  the  probability  density  function  implied  by  Equation  1,  it  is  necessary  to  rewrite 
the  left  sidP  of  the  equation  as 


N  •  t2 

N  ■  E  (failures  by  t2 ) 


K  t: 


1/3 


(2) 


where 


E(failuris  t>y  t3 )  is  the  expected  number  of  failures  from  time  zero  to  time  t2 


The  x ported  number  of  failures  is  actually  the  number  of  repairs  for  a  renewable 
proris.-.  Therefore,  equations  for  the  reliability  of  renewable  units  can  be  used  to  derive  a 
density  function  from  Equation  7. 

For  a  renewable  process,  the  cumulative  number  of  renewals  (or  repairs  or  failures)  is 
given  by  the  “renewal  function”  II(t).  From  Equation  2,  the  renewal  function  is  given  by 

ll(tj)  -  E(failures  by  t3)  -  t22/3/K 


B~1 


(3) 


It  can  be  shown*  that  the  renewal  function  is  bounded  by 


1  -  U(t)  <  H(t) 


l-K(t) 

R(t) 


where 


R(t)  -  l-F(t)  is  the  reliability  function 
If  we  use  the  relation 


1  -  R(t)  <  -  fin  R(t)  < 


l-R(t) 

R(t) 


(4) 


(5) 


for  0  <  R(t)  <  1,  then 


H(ta)  *»  -  In  R(tj ) 


(6) 


provided  the  outer  terms  in  inequalities  4  and  5  are  nearly  identical.  For  highly  reliable 
systems,  R(t2 )  will  be  nearly  equal  to  one.  Therefore, 


1  -  R(t,) 


1  -R(tj) 
R(tj) 


(7) 


Thus  approximation  6  is  acceptable  for  this  analysis.** 

Combining  Equations  3  and  6,  we  obtain 

8n  R(tj)  -  -tj2/3/K  (8) 


and 


R(t2)  - 


P-t:2/3/r 


(9) 


Equation  9  can  be  written  in  the  forn:  of  the  reliability  function  for  the  Wcibull 
distribution: 


--a  t,2/3 
R\v(li)  “  e  2 
2 

with  n  “  —  (19) 

3K 


*B.V.  Gnrdo;il:o,  Yu.  K.  Kclycycr,  pml  A.D.  Solovyev,  Krtha.inticrt  l'clho:'s  of  Reliability  Theory , 
Academic  IYoks,  New  York,  1900,  pp  93-102. 

**For  cxcir-;,!: ,  if  \vc  lire  Equation  10  viih  t2  “  14,000  Iiox»r s,  tiicn  l"R(f2)  **  0.007758  rncl 
ll-R(t2)l/R(’l2)-  0.007819. 


B~2 


The  resulting  density  function  and  cumulative  distribution  are  then,  respectively. 


Mt,)  -  „„ 

and 

Fw(t2 )  -=  i  _  <r  3/2  at2  2/3  (12) 

2 

Thus  the  regression  analysis  yields  a  shape  parameter  of  2/3  and  a  scale  parameter  of  - — 

3K 

When  Equation  10  is  rewritten  in  terms  of  the  derived  coefficients,  the  reliability  function 
and  the  range  of  parameter  values  supported  by  the  data  become 

R„(t.)  ■  .-‘’2/3'k 


where 


K  =  7G  877  ‘  e^O.025  Sc  +  O.OOOS5  +  0.0034  tj 
Sc  »  screening  score,  0  <  Sc  <  70 
Sa  -  sampling  score,  0  <  Sa  <  2,438 
t!  -  system  burn-in  lime,  0  <  t(  <311  hours 
t2  -  Geld  operate  time,  9G  <  t2  <  14,600  hours 
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